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1 Introduction 

The early Universe was hot and dense behaving as a 
cosmic nuclear reactor during the first twenty minutes 
of its evolution. It was, however, a "defective" nuclear 
reactor, expanding and cooling very rapidly. As a re- 
sult, only a handful of the lightest nuclides were synthe- 
sized before the density and temperature dropped too 
low for the nuclear reaction rates to compete with the 
universal expansion rate. After hydrogen (^H = pro- 
tons) the next most abundant element to emerge from 
the Big Bang is helium ^He = alpha particles). Iso- 
topes of these nuclides (deuterium and helium-3) are 
the next most abundant primordially. Then there is 
a large gap to the much lower abundance of lithium- 
7. The relative abundances of all other primordially- 
produced nuclei are very low, much smaller than their 
locally observed (or, presently observable!) abundances. 
After a brief description of the early evolution of the 
Universe emphasizing those aspects most relevant to 
primordial, or "Big Bang" nucleosynthesis (BBN), the 
predicted abundances of the light nuclides will be pre- 
sented as a function of the one "free" parameter (in the 
simplest, "standard" model: SBBN), the nucleon (or 
"baryon" ) abundance. Then, each element will be con- 
sidered in turn in a confrontation between the predic- 
tions of SBBN and the observational data. At present 
(Summer 1999) there is remarkable agreement between 
the SBBN predictions of the abundances of four nuclides 
(D, "^He, *IIe, and ^Li) and their primordial abundances 
inferred from the observations. However, there are some 
hints that this concordance of the hot big bang model 
may be imperfect, so we will also explore some varia- 
tions on the theme of the standard model with regard 
to their modifications of the predicted primordial abun- 
dances of the light elements. 

In the simplest, standard, hot big bang model the 
currently observed large-scale isotropy and homogene- 
ity of the Universe is assumed to apply during ear- 
lier epochs in its evolution. Given the currently ob- 
served universal expansion and the matter and radia- 
tion (CBR: "cosmic background radiation", the 2.7K 
"black body radiation") content, it is a straightforward 
application of classical physics to extrapolate back to 
earlier epochs in the history of the Universe. At a time 
of order 0.1 s after the expansion began the Universe 
was filled with a hot, dense plasma of particles. The 
most abundant were photons, electron-positron pairs, 
particle-antiparticle pairs of all known "flavors" of neu- 
trinos [i/e, v^, and i/r), and trace amounts of neutrons 



and protons ("nucleons" or "baryons"). At such early 
times the thermal energy of these particles was very 
high, of order a few MeV. With the exception of the 
nucleons, it is known or assumed that all the other par- 
ticles present were extremely relativistic at this time. 
Given their high energies (and velocities close to, or ex- 
actly equal to the speed of light) and high densities, 
the electroweak interactions among these particles were 
sufficiently rapid to have established thermal equilib- 
rium. As a result, the numbers and distributions (of 
momentum and energy) of all these particles is accu- 
rately predicted by well-known physics. 

2 Nucleosynthesis In The Early 
Universe 

The primordial yields of light elements are determined 
by the competition between the expansion rate of the 
Universe (the Hubble parameter H) and the rates of the 
weak and nuclear reactions. It is the weak interaction, 
interconverting neutrons and protons, that largely de- 
termines the amount of ^Hc which may be synthesized, 
while detailed nuclear reaction rates regulate the pro- 
duction (and destruction) of the other light elements. 
In the standard model of cosmology the early expan- 
sion rate is fixed by the total energy density p, 



H^ = 8ttGp/3, 



(1) 



where G is Newton's gravitational constant. In the 
standard model of particle physics the early energy den- 
sity is dominated by the lightest, relativistic particles. 
For the epoch when the Universe is a few tenths of a 
second old and older, and the temperature is less than 
a few MeV, 

/9 = p^ + Pe + N^/3^, (2) 

where p-y, pe, and pi, are the energy densities in photons, 
electrons and positrons, and massless neutrinos and an- 
tineutrinos (one species), respectively; A^^ is the num- 
ber of massless (or, very light: m^ <C 1 MeV) neutrino 
species which, in standard BBN, is exactly 3. In con- 
sidering variations on the theme of the standard model, 
it is useful to allow N^, to differ from 3 to account for 
the presence of "new" particles and/or any suppression 
of the standard particles {e.g., if the t neutrino should 
have a large mass). Since the energy density in rela- 
tivistic particles scales as the fourth power of the tem- 
perature, the early expansion rate scales as the square 
of the temperature with a coefBcient that depends on 



the number of different relativistic species. The more 
such species, the faster the Universe expands, the ear- 
her (higher temperature) will the weak and nuclear re- 
actions drop out of equilibrium. It is useful to write the 
total energy density in terms of the photon energy den- 
sity and g, the equivalent number of relativistic degrees 
of freedom {i.e., helicity states, modulo the different 
contributions to the energy density from fermions and 
bosons), 

p ^ {9/2)p^. (3) 

In the standard model at T -- 1 MeV, gsM = 43/4. 
Account may be taken of additional degrees of freedom 
by comparing their contribution to p to that of one ad- 
ditional light neutrino species. 



Ap = ptot - PSM = ^N^p^. 



(4) 



If the early energy density deviates from that of the 
standard model, the early expansion rate (or, equiva- 
lently, the age at a fixed temperature) will change as 
well. The "speed-up" factor ^ = H/Hsm may be re- 
lated to AN^ by. 



e=(p/P5M)^/' = (l + 7AiV./43)i/2 



(5) 



As we'll see shortly, the ''He abundance is very sensi- 
tive to the early expansion rate while the abundances of 
the other light nuclides depend mainly on the nuclear 
reaction rates which scale with the nucleon (baryon) 
density. Since the baryon density is always changing 
as the Universe expands, it is convenient to distinguish 
between models with different baryon densities using a 
dimensionless parameter which is either conserved or, 
changes in a known and calculable fashion. From the 
very early Universe till now the number of baryons in 
a comoving volume has been preserved and the same is 
roughly true for photons since the end of BBN. There- 
fore, the ratio of number densities of baryons (riB) and 
photons (n-y) provides just such a measure of the uni- 
versal baryon abundance. 



jj = (nB/%)o 



Vio = 10^°V- 



(6) 



The universal density of photons at present (through- 
out this article the present epoch is indicated by the 
subscript '0') is dominated by those in the CBR (for 
To = 2.73 K, n-yo = 412 cm~'^) so that the baryon den- 
sity parameter J7b = (pb/Pc)o, the ratio of the present 
baryon density (pe) to the present critical density (pc), 
may be related to rj and the present value of the Hubble 
parameter Hq = lOO/i kms~^Mpc~^, 



ryio = 273f^Bft- • 



(7) 



It should be noted that prior to electron-positron an- 
nihilation there were fewer photons in every comoving 
volume (by a factor very close to 4/11); this is auto- 
matically accounted for in all numerical BBN codes. It 
is simply a matter of consensus and convenience that 
the baryon abundance is quoted in terms of its present 
value. 

In SBBN {i.e., N^ = 3) the abundances of the light 
nuclides synthesized primordially depend on only one 
"free" parameter, rj. SBBN is thus "overconstrained" 
since one value (or, a narrow range of values set by the 
observational and theoretical reaction rate uncertain- 
ties) of rj must account consistently for the primordial 
abundances of D, ■^He, ''He and ^Li. At the same time 
this value/range of rj must be consistent with current 
estimates of (or, bounds to) the present baryon density. 
For these reasons BBN is one of the key pillars support- 
ing the edifice of the standard model of cosmology and, 
it is the only one which offers a glimpse of the earliest 
evolution of the Universe. In the following we'll first 
identify the key landmarks in the first 20 minutes in 
the evolution of the Universe in order to identify the 
physical processes responsible for determining the pri- 
mordial abundances of the light nuclides. Then, after 
presenting the SBBN predictions (as a function of tj; see 
Fig. 1) we will review the current status of the observa- 
tional data, as well as the steps necessary in order to go 
from "here and now" to "there and then" when using 
the data to infer the true primordial abundances. Then 
we will be in a position to assess the consistency of the 
standard model. 

2.1 Weak Equilibrium and the ^He 
Abundance 

Consider now those early epochs when the Universe was 
only a few tenths of a second old and the radiation 
filling it was at a temperature (thermal energy) of a 
few MeV. According to the standard model, at those 
early times the Universe was a hot, dense "soup" of 
relativistic particles (photons, e^ pairs, 3 "flavors" (e, 
fi, t) of neutrino-antineutrino pairs) along with a trace 
amount (at the level of a few parts in 10""') of neutrons 
and protons. At such high temperatures and densities 
both the weak and nuclear reaction rates are sufficiently 
rapid (compared to the early Universe expansion rate) 
that all particles have come to equilibrium. A key con- 
sequence of equilibrium is that the earlier history of 
the evolution of the Universe is irrelevant for an under- 
standing of BBN. When the temperature drops below a 
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Figure 1: The predicted primordial abundances as a 
function of r/. Yp is the ^He mass fraction while j/2P, 
2/3Pj 2/7P are the number density ratios to hydrogen of 
D, "^He, and *He respectively. 



few McV the weakly interacting neutrinos effectively de- 
couple from the photons and c pairs, but they still play 
an important role in regulating the ncutron-to-proton 
ratio. 

At high temperatures, neutrons and protons are 
continuously intcrconvcrting via the weak interactions: 
n + e+^^p + j^e, n + j/g^^P + e^, and n ^^ p + e^ + i>e- 
When the interconversion rate is faster than the expan- 
sion rate, the ncutron-to-proton ratio tracks its equilib- 
rium value, decreasing exponentially with temperature 
{n/p = e~^™/^, where Am = 1.29 MeV is the neutron- 
proton mass difference). A comparison of the weak rates 
with the universal expansion rate reveals that equilib- 



rium may be maintained until the temperature drops 
below ^ 0.8 MeV. When the interconversion rate be- 
comes less than the expansion rate, the n/p ratio ef- 
fectively "freezes-out" (at a value of sa 1/6), thereafter 
decreasing slowly, mainly due to free neutron decay. 

Although n/p freeze-out occurs at a temperature 
below the deuterium binding energy, Eb = 2.2 MeV, 
the first link in the nucleosynthetic chain, p + n ^ 
D + 7, is ineffective in jump-starting BBN since the 
photodcstruction rate of deuterium (oc n-^e^^'^'^) is 
much larger than the deuterium production rate (ex ub) 
due to the very large universal photon-to-baryon ratio 
(> 10^). Thus, the Universe must "wait" until there are 
so few sufficiently energetic photons that deuterium be- 
comes effectively stable against photodissociation. This 
occurs for temperatures ^ 80 keV, at which time neu- 
trons are rapidly incorporated into ^He with an effi- 
ciency of 99.99%. This efficiency is driven by the tight 
binding of the ''He nucleus, along with the roadblock to 
further nucleosynthesis imposed by the absence of a sta- 
ble nucleus at mass-5. By this time (T ^ 80 keV), the 
n/p ratio has dropped to ~ 1/7, and simple counting 
(2 neutrons in every ^He nucleus) yields an estimated 
primordial *He mass fraction 



Yp 



2{n/p) 
[1 + (n/p)] 



0.25. 



(8) 



As a result of its large binding energy and the gap at 
mass-5, the primordial abundance of ^He is relatively 
insensitive to the nuclear reaction rates and, therefore, 
to the baryon abundance (rj) . As may be seen in Figure 
1, while 1] varies by orders of magnitude, the predicted 
"'He mass fraction, Yp , changes by factors of only a few. 
Indeed, for 1 < ??io < 10, 0.22 < Yp < 0.25. As may be 
seen in Figures 1 and 2, there is a very slight increase in 
Yp with rj. This is mainly due to BBN beginning earlier, 
when there arc more neutrons available to form ^He, if 
the baryon-to-photon ratio is higher. The increase in 
Yp with rj is logarithmic; over most of the interesting 
range in rj, AYp « O.OIA77/77. 

The ''He abundance is, however, sensitive to the 
competition between the universal expansion rate (H) 
and the weak interaction rate (intcrconvcrting neutrons 
and protons) . If the early Universe should expand faster 
than predicted for the standard model, the weak inter- 
actions will drop out of equilibrium earlier, at a higher 
temperature, when the n/p ratio is higher. In this 
case, more neutrons will be available to be incorpo- 
rated into ''He and Yp will increase. Numerical cal- 
culations show that for a modest speed-up {AN^^ 1), 
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Figure 2: The predicted "^He abundance (solid curve) 
and the 2a theoretical uncertainty 0. The horizon- 
tal lines show the range indicated by the observational 
data. 

AYp w 0.013AA''i,. Hence, constraints on Yp (and rf) 
lead directly to bounds on AA^i, and, on particle physics 
beyond the standard model ||l| . 

It should be noted that the uncertainty in the BBN- 
predicted mass fraction of ''He is very small and almost 
entirely dominated by the (small) uncertainty in the 
n — p interconversion rates. These rates may be "nor- 
malized" through the neutron lifetime, r„, whose cur- 
rent standard value is 887±2 s (actually, 886.7±1.9 s). 
To very good accuracy, a 1 s uncertainty in t„ corre- 
sponds to an uncertainty in Yp of order 2 x 10^''. At 
this tiny level of uncertainty it is important to include 
finite mass, zero- and finite-temperature radiative cor- 
rections, and Coulomb corrections to the weak rates. 
However, within the last few years it emerged that the 
largest error in the BBN-prediction of Yp was due to 
a too large time-step in the numerical code. With this 
now under control, it is estimated that the residual the- 
oretical uncertainty (in addition to that from the un- 
certainty in T„) is of the order of 2 parts in 10'*. Indeed, 
a comparison of two major, independent BBN codes re- 
veals agreement in the predicted values of Yp to 0.0001 
± 0.0001 over the entire range 1 < tjiq < 10. In Figure 2 
is shown the BBN-predicted ''He mass fraction, Yp , as 
a function of 77; the thickness of the band is the ±2cr 
theoretical uncertainty. For tjiq > 2 the la theoretical 
uncertainty in Yp is ^ 6 x 10~*. As we will soon see, 
the current observational uncertainties in Yp are much 
larger (see, also. Fig. 2). 



2.2 Deuterium - The Ideal Baryometer 

As may be seen in Fig. 1, the deuterium abundance 
(the ratio, by number, of deuterium to hydrogen: here- 
inafter, (D/H)p = j/2p) is a monotonic, rapidly decreas- 
ing function of the baryon abundance 77. The reason for 
this behavior is easily understood. Once BBN begins in 
earnest, when the temperature drops below '^ 80 keV, 
D is rapidly burned to '^H, '^He and *He. The higher the 
baryon abundance, the faster the burning and the less D 
survives. For 7710 in the "interesting" range 1 - 10, y2P 
decreases with the ~ 1.6 power of 77. As a result, a 10% 
error in j/2P corresponds to only a 6% error in 77. This 
strong dependence of y2P on r^io, combined with the 
simplicity of the evolution of D/H in the epochs follow- 
ing BBN, is responsible for the unique role of deuterium 
as a baryometer M . Because almost all the relevant re- 
action cross sections are measured in the laboratory at 
energies comparable to those of BBN, the theoretical 
uncertainties in the BBN-predicted abundance of deu- 
terium is quite small, 8 - 10% for most of the interesting 
rj range shown in Fig. 3. 
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Figure 3: The predicted D/H abundance (solid curve) 
and the 2a theoretical uncertainty g]. The horizontal 
lines show the range indicated by the observational data 
for both the high D/H (upper two lines ) and low D/H 
(lower two lines). 

Deuterium and helium-4 are complementary, form- 
ing the crucial link in testing the consistency of BBN in 
the standard model. While the primordial-D abundance 
is very sensitive to the baryon density, the primordial- 
*He abundance is relatively insensitive to rj. Deuterium 
provides a bound on the universal baryon density while 



helium-4 constrains the early expansion rate of the Uni- 
verse, ofFerring bounds on particle physics beyond the 
standard model. 

2.3 Helium-3 — Complicated Evolution 

As may be seen in Fig. 1, the predicted primordial abun- 
dance of ^He behaves similarly to that of D, decreasing 
monotonically with rj. Again, the reason is the same: 
^He is being burned to the more tightly bound nucleus 
■^He. The higher the nucleon abundance, the faster the 
burning and the less ^He survives. In contrast to the 
behavior of D/H versus ij, the decrease of ^He/H with 
increasing rj is much slower. This is simply a reflection 
of the much tighter binding of ^He compared to D. Al- 
though the BBN predictions for the abundance of ^He 
have similarly small uncertainties (8 - 10%) to those of 
D, it is much more difficult to exploit the observations 
of '^He to test and constrain BBN. The complication 
is the evolutionary history of the '^He abundance since 
BBN. 

Although any deuterium cycled through stars is 
burned to ^He during the stars' pre-main sequence evo- 
lution, '^He will survive in the cooler stellar exteriors 
while being destroyed in the hotter interiors. For the 
more abundant lower mass stars which are cooler, a 
larger fraction of prestellar ^He (along with the ^He 
produced from prestellar D) survives. Indeed, for suf- 
ficiently low-mass stars (less than a few solar masses) 
incomplete burning actually leads to a buildup of newly 
synthesized ^He (to be contrasted with the prestellar D 
and '^He) which may - or may not - be returned to 
the interstellar medium. In fact, some planetary neb- 
ulae are observed to be highly enriched in ^He. So, 
the evolution of ■^He is complex with stellar destruc- 
tion competing with primordial and stellar production. 
Indeed, if all low mass stars were as prolific produc- 
ers of '^He as indicated by some planetary nebulae, the 
solar system and local interstellar medium abundances 
of ^He should far exceed those inferred from observa- 
tions. Thus, at least some low mass stars must be net 
destroyers of '^He. Given this necessarily complex and 
uncertain picture of production, destruction and sur- 
vival, it is difficult to use current observational data to 
infer the primordial abundance of ■^He. Unless and until 
^He is observed in high rcdshift {i.e., early Universe), 
low metallicity {i.e., nearly unevolved) systems, it will 
provide only a weak check on the consistency of BBN. 



2.4 Lithium-7 The Lithium Valley 

The trend of the BBN-predicted primordial abundance 
of lithium (almost entirely ^Li) with r/ is more 'interest- 
ing' than that of the other light nuclides (see Figs. 1 & 
4). The "lithium valley", centered near 7710 w 2 — 3, is 
the result of the competition between production and 
destruction in the two paths to mass-7 synthesis. At rel- 
atively low baryon abundance (7710^ 2) mass-7 is mainly 
synthesized as ^Li via '^H +'*IIe — >^Li + 7. As the 
baryon abundance increases at low 77, ^Li is destroyed 
rapidly by {p,^)2a reactions. Hence the decrease in 
(Li/II)p = yjp with increasing rj seen (at low 77) in 
Figs. 1 & 4. Were this the only route to primordial 
synthesis of mass-7, this monotonic trend would con- 
tinue, similar to those for D and "^He. However, mass-7 
may also be synthesized via "^He -|-''^He ^-^Be -f- 7. The 
^Be will later capture an electron to become ^Li. This 
channel is very important because it is much easier to 
destroy ^Li than ^Be. As a result, for relatively high 
baryon abundance (7/10 >; 3) this latter channel domi- 
nates mass-7 production and yjp increases with increas- 
ing 77. 
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Figure 4: The predicted ^Li abundance (solid curve) 
and the 2cr theoretical uncertainty M. The horizon- 
tal lines show the range indicated by the observational 
data. 

As may been seen in Figure 4, the BBN-predicted 
uncertainties for 7/7P are much larger than those for D, 
■^He, or ^He. In the interval 1 < 7710 < 10 the la uncer- 
tainties are typically ~ 20%, although in a narrow range 
of 7] near the bottom of the 'valley' they are somewhat 
smaller (- 12 - 15%). 



3 From Here And Now To There 
And Then 

To test the consistency of SBBN requires that we con- 
front the predictions with the primordial abundances of 
the Hght nucHdes which, however, are not observedhut, 
rather, are inferred from observations. The path from 
the observational data to the primordial abundances 
is long and twisted and often fraught with peril. In 
addition to the usual statistical and systematic uncer- 
tainties, it is crucial to forge a connection from "here 
and now" to "there and then"; i.e., to relate the de- 
rived abundances to their primordial values. It is indeed 
fortunate that each of the key elements is observed in 
different astrophysical sites using very different astro- 
nomical techniques. Also, the corrections for chemical 
evolution differ among them and, even more important, 
they can be minimized. For example, deuterium (and 
hydrogen) is mainly observed in cool, neutral gas (so 
called Hi regions) via UV absorption from the atomic 
ground state (the Lyman series), while radio telescopes 
allow helium-3 to be studied utilizing the analog of the 
hydrogen 21 cm line for singly-ionized '^He in regions 
of hot, ionized gas (so called Hii regions). Helium-4 is 
probed using the emission from its optical recombina- 
tion lines formed in Hii regions. In contrast, lithium is 
observed in the absorption spectra of warm, low-mass 
halo stars. With such different sites, with the mix of 
absorption/emission, and with the variety of telescopes 
and different detectors involved, the possibility of corre- 
lated errors biasing the comparison with the predictions 
of BBN is unlikely. This favorable situation extends to 
the obligatory evolutionary corrections. For example, 
although until recently observations of deuterium were 
limited to the solar system and the Galaxy, mandat- 
ing uncertain evolutionary corrections to infer the pre- 
galactic abundance, the Keck and Hubble Space tele- 
scopes have begun to open the window to deuterium in 
high-redshift, low-metallicity, nearly primordial regions 
(Lyman-a clouds). Observations of *He in chemically 
uncvolved, low-metallicity {^ 1/50 of solar) extragalac- 
tic H II regions permit the evolutionary correction to 
be reduced to the level of the statistical uncertainties. 
The abundances of lithium inferred from observations 
of the very metal-poor halo stars (one-thousandth of 
solar abundance and even lower) require almost no cor- 
rection for chemical evolution. On the other hand, as 
noted earlier, the status of helium-3 is in contrast to 
that of the other light elements. For this reason, '^He 



will not be used quantitatively in this article. 

The currently very favorable observational and evo- 
lutionary situation for the nuclides produced during 
BBN is counterbalanced by the likely presence of sys- 
tematic errors in the path from observations to primor- 
dial abundances. By their very nature, such errors are 
difhcult - if not impossible - to quantify. In the key 
case of deuterium there is a very limited set of the most 
useful data. As a result, and although cosmological 
abundance determinations have taken their place in the 
current "precision" era of cosmology, it is far from clear 
that the present abundance determinations are truly 
"accurate". Thus, the usual caveat emptor applies to 
any conclusions drawn from the subsequent compari- 
son between the predictions and the data. With this 
caution in mind the current status of the data will be 
surveyed in order to infer "reasonable" ranges for the 
primordial abundances of the key light elements. 

3.1 Deuterium 

Since deuterium is completely burned away whenever 
it is cycled through stars, and there are no astrophys- 
ical sites capable of producing deuterium in anywhere 
near its observed abundance 0, any D-abundance de- 
rived from observational data provides a lower bound 
to its primordial abundance. Thus, without having to 
correct for Galactic evolution, the 'here-and-now' deu- 
terium abundance inferred from UV observations along 
12 lines-of-sight in the local interstellar medium (LISM) 
bounds the 'there-and-then' primordial abundance from 
below ((D/H)p >(D/H)lism = (1.5 ±0.1) x 10"^) ^s 
may be seen from Figs. 1 and 2, any lower bound to pri- 
mordial D will provide an upper bound to the baryon- 
to-photon ratio M. 

Solar system observations of '^He permit an inde- 
pendent, albeit indirect determination of the pre-solar 
system deuterium abundance p] . This estimate of the 
Galactic abundance some 4.5 Gyr ago, while somewhat 
higher than the LISM value, has a larger uncertainty 
(D/H = (2.1 ± 0.5) X 10-5 H). Within the uncer- 
tainty it is consistent with the LISM value suggesting 
there has been only modest destruction of deuterium 
in the last 4.5 Gyr. There is also a recent measure- 
ment of deuterium in the atmosphere of Jupiter us- 
ing the Galileo Probe Mass Spectrometer |Q, (D/H 
= (2.6 ±0.7) X 10-5). 

To further exploit the solar system and/or LISM 
deuterium determinations to constrain/estimate the 
primordial abundance would require corrections for the 



Galactic evolution of D. Although the simplicity of the 
evolution of deuterium (it is only destroyed) suggests 
that such correction might be very nearly indepen- 
dent of the details of specific chemical evolution mod- 
els, large differences remain between different estimates. 
It is therefore fortunate that data on D/H in high- 
redshift (nearly primordial), low-metallicity (nearly un- 
evolved) Lyman-a absorbers has become available in 
recent years (§,§)■ It is expected that such systems 
still retain their original, primordial deuterium, undi- 
luted by the deuterium-depicted debris of any signif- 
icant stellar evolution. That's the good news. The 
bad news is that, at present, D-abundance determi- 
nations are claimed for only three such systems, and 
that the abundances inferred for two of them (along 
with a limit to the abundance for a third) appear to 
be inconsistent with the abundance estimated for the 
remaining one. Here we have a prime illustration of 
"precise", but possibly inaccurate cosmological data. 
Indeed, there is a serious obstacle inherent to using ab- 
sorption spectra to measure the deuterium abundance 
since the isotope-shifted deuterium absorption lines are 
indistinguishable from the corresponding lines in suit- 
ably velocity-shifted (—81 kms~^) hydrogen. Such "in- 
terlopers" may have been responsible for some of the 
early claims [|| of a "high" deuterium abundance |]lO| . 
Indeed, an interloper may be responsible for the one 
surviving high-D claim ||ll|. At present it seems that 
only three good candidates for nearly primordial deu- 
terium have emerged from ground- and space-based 
observations. It may be premature to constrain cos- 
mology on the basis of such sparse data. Nonethe- 
less, the two "low-D" systems suggest a primordial deu- 
terium abundance consistent with estimates of the pre- 
Galactic value inferred from LISM and solar system 
data ((D/H)p = (3.4 ±0.25) x IQ-^ ||). To illustrate 
the confrontation of cosmological theory with observa- 
tional data, this D-abundance will be adopted in the 
following. However, the consequences of choosing the 
"high-D" abundance ((D/H)p = (20 ± 5) x 10-^ 0) 
will also be discussed. 

3.2 Helium-4 

As the second most abundant nuclide in the Universe 
(after hydrogen), the abundance of "'He can be deter- 
mined to high accuracy at sites throughout the Uni- 
verse. However, as stars evolve they burn hydrogen to 
helium and the ^He in the debris of stellar evolution 
contaminates any primordial ^He. Since any attempt 



to correct for stellar evolution will be inherently uncer- 
tain, it is sensible to concentrate on the data from low- 
metallicity sites. Extragalactic regions of hot, ionized 
gas (Hii regions) provide such sites, where the helium is 
revealed via emission lines formed when singly and dou- 
bly ionized helium recombines. As with deuterium, cur- 
rent data provide ambiguous estimates of the primordial 
helium abundance. Since the differences (AY = 0.010) 
are larger than the statistical uncertainties (< ±0.003), 
systematic errors likely dominate. Among the currently 
most likely sources of such errors are uncertain cor- 
rections for coUisional excitation in helium, uncertain 
corrections for unseen neutral helium and/or hydrogen, 
and underlying stellar absorption (leading to an under- 
estimate of the true strength of the helium emission 
lines). In contrast, since the most metal-poor of the ob- 
served regions have metallicities of order 1/50 - 1/30 of 
solar, the extrapolation from the lowest metallicity re- 
gions to truly primordial introduces an uncertainty no 
larger than the statistical error. 

Using published data ( |1J, ^ ) for 40 low-metallicity 
regions, one group |]l| finds: Yp = 0.234 ± 0.003. 
In contrast, from an independent data set of 45 low- 
metallicity regions another group |lq| infers Yp = 
0.244 ± 0.002. Clearly, these results are statistically 
inconsistent. It is crucial that high priority be assigned 
to further Hii region observations to estimate/avoid the 
systematic errors. Until then, since the error budget for 
Yp is likely dominated by systematic rather than sta- 
tistical uncertainties, in what follows, a generous range 
for Yp wiU be adopted: 0.228 < Yp < 0.248. 

3.3 Lithium-7 

Cosmologically interesting lithium is observed in the 
spectra of nearly 100 very metal-poor halo stars ll^ ESJ . 
These stars are so metal-poor they provide a sample 
of more nearly primordial material than anything ob- 
served anywhere else in the Universe; the most metal- 
poor among them have heavy element abundances less 
than one-thousandth of the solar metallicity. However, 
these halo stars are also the oldest objects in the Galaxy 
and, as such, have had the most time to modify their 
surface abundances. So, even though any correction for 
Galactic evolution modifying their lithium abundances 
may be smaller than the statistical uncertainties of a 
given measurement, the systematic uncertainty asso- 
ciated with the dilution and/or destruction of surface 
lithium in these very old stars could dominate the er- 
ror budget. There could be additional errors associ- 



ated with the modeling of the surface layers of these 
cool, low-metallicity, low-mass stars needed to derive 
abundances from absoprtion-line spectra. It is also 
possible that some of the Li observed in these stars 
is non- primordial {e.g., that some of the observed Li 
may have been produced post-BBN by spallation reac- 
tions (the breakup of C, N, and O nuclei into nuclei of 
Li, Be, and B) or fusion reactions (a -I- a to form ^Li 
and ^Li) in cosmic-ray collisions with gas in the ISM). 
In a recent analysis (20|, it is argued that as much as 
^ 0.2 dex of the observed lithium abundance, A(Li) 
= 12-flog(Li/H), could be post-primordial in origin. 

The very large data set of lithium abundances mea- 
sured in the warmer (T > SSOOiiT), more metal-poor 
([Fe/H] < —1.3) halo stars define a plateau (the "Spite- 
plateau" |1^) in the lithium abundance - metallicity 
plane. Depending on the choice of stellar-temperature 
scale and stellar atmosphere model, the abundance level 
of the plateau is: A(Li) = 2.2 ± 0.1, with very little dis- 
persion in abundances around this plateau value. The 
small dispersion provides an important constraint on 
models which attempt to connect the present surface 
lithium abundances in these stars to the original lithium 
abundance in the gas out of which these stars were 
formed some 10 - 15 Gyr ago. "Standard" (i.e., non- 
rotating) stellar models predict almost no lithium de- 
pletion and, therefore, are consistent with no disper- 
sion about the Spite-plateau. Although early work on 
mixing in models of rotating stars was very uncertain, 
recently stellar models have been constructed which re- 
produce the angular momentum evolution observed for 
the much younger, low-mass open cluster stars. These 
models have been applied to the study of lithium deple- 
tion in main sequence halo stars. A well-defined lithium 
plateau with modest scatter and a small population 
of "outliers" (overdepleted stars), consistent with the 
current data, is predicted for depletion factors between 



0.2 dex and 0.4 dex flQ 



To err on the side of caution, a generous range for 
the plateau abundance, 2.1 < A(Li) < 2.3, is adopted. 
If depletion is absent, this range is consistent with the 
primordial lithium "valley" (see Fig. 4). For depletion 
> 0.2 dex, the consistent primordial lithium abundances 
bifurcate and move up into the "foothills", although 
a non-negligible contribution from post-BBN lithium 
could move the primordial abundance back down again 
(Fig. 4). 



4 Confrontation Of Theory With 
Data 

In the context of the "standard" model (three families of 
massless, or light, two-component neutrinos), the pre- 
dictions of BBN (SBBN) depend on only one free pa- 
rameter, the nucleon-to-photon ratio rj. The key test 
of the standard, hot, big bang cosmology is to assess if 
there exists a unique value or range of rj for which the 
predictions of the primordial abundances are consistent 
with the light element abundances inferred from the 
observational data. From a statistical point of view it 
might be preferrable to perform a simultaneous fit of the 
inferred primordial abundances of D, ■^He, "^He, and ''Li 
to the SBBN predictions. In this manner the "best fit" 
77, along with its probability distribution may be found, 
and the "goodness-of-fit" assessed [|l|. However, since 
systematic uncertainties most likely dominate observa- 
tional errors at present, the value of this approach is 
compromised. An alternate approach is adopted here. 

As emphasized earlier, deuterium is an ideal bary- 
ometer. As a first step the primordial abundance of 
deuterium inferred from observations at high redshift 
will be compared with the SBBN prediction to iden- 
tify a consistent range for rj. Then, given this range, 
the SBBN abundances of ''He and ^Li are predicted 
and these are compared to the corresponding primor- 
dial abundances derived from the observational data. 
The challenge is to see if the D-identified range for 77 
leads to consistent predictions for ''He and ^Li. Recall 
that due to its complicated evolutionary history, it is 
difficult to use ^He to test and constrain SBBN. Fur- 
thermore, another consistency test is to compare the 
SBBN-inferred rj range with the present baryon den- 
sity derived from non-BBN observations and theory. Is 
our model for the very early evolution of the Universe 
consistent with the present Universe? 

From the two well observed, high redshift absorption 
line systems with "low-D" , the estimate adopted for the 
primordial-D abundance is: (D/H) p = 2.9 — 4.0 x 10~^ 
(see Fig. 3). Also shown for comparison in Figure 3 
is the allowed range of the primordial deuterium abun- 
dance suggested by the "high-D" abundance inferred 
from observations of one lower redshift absorption-line 
system. With allowance for the ^ 8% uncertainty in the 
theoretically predicted abundance, the favored range 
(low-D) for T] is quite narrow: 7710 = 5.1 ± 0.36. It 
is clear from Figure 4 that for the baryon abundance 
in this range, the BBN-predicted lithium abundance is 



entirely consistent with the Spite-plateau value, even if 
the plateau were raised by ~ 0.2 dex to allow for mod- 
est stellar destruction/dilution or lowered by a similar 
amount due to post-BBN production. For this narrow 
range in rj the predicted ^He mass fraction varies very 
httle. For ryio « 5, AYp « O.OlOAry/Ty, so that in- 
cluding the error in the predicted abundance, Yp = 
0.247 ± 0.001. As may be verified from Figure 2, this is 
within (albeit at the high end of) the range allowed by 
the data from the low metallicity, extragalactic H ii re- 
gions. Given the current uncertainties in the primordial 
abundances, SBBN is consistent with "low-D/high-?7" . 

The significance of this concordance cannot be un- 
derestimated. A glance at Fig. 1 provides a reminder of 
the enormous ranges for the predicted primordial abun- 
dances. That the simplest hot, big bang cosmological 
model can account for ("predict" !) 3 independent abun- 
dances (4 with ■^He; although ^He hasn't been employed 
in this comparison, its predicted abundance is consis- 
tent with extant observational data) by adjusting only 
one free parameter {r]) is a striking success. The the- 
ory, which is in principle falsifiable, has passed the test. 
It needn't have. Indeed, future observational data cou- 
pled to better understanding of systematic errors may 
provide new challenges. For example, if in the future it 
should be determined that the primordial helium mass 
fraction were lower than Yp = 0.245, this would be in- 
consistent (within the errors) with the "low-D/high-77" 
range derived above. Similarly, if the best estimate for 
the D-determined 77 range changed, the comparison per- 
formed above should be repeated. With this in mind, 
what of the "high-D/low-Ty" range which has been set 
aside in the current comparison? 

If, in contrast to the deuterium abundance adopted 
above, the true value were higher, (D/II)p = 10 — 30 x 
10~^, the SBBN- favored range in rj would be lower (see 
Fig. 2). Accounting for the ^ 8% uncertainty in the 
theoretically predicted abundance, ryio = 1.7 ± 0.28. 
Inspection of Figures 2-4 reveals that as along as 
''710 ^ 1-1 — 1-3, consistency with ^He and ''Li can be 
obtained. Hence, for "high-D" as well, the standard 
model passes the key cosmological test. 

4.1 Comparison Of BBN With Non- 
BBN Baryon Density Estimates 

Having established the internal consistency of primor- 
dial nucleosynthesis in the standard model, it is neces- 
sary to proceed to the next key test. Does the nucleon 
abundance inferred from processes which occurred dur- 



ing the first thousand seconds of the evolution of the 
Universe agree with estimates/bounds to the nucleon 
density in the present Universe? 

It is a daunting task to attempt to inventory the 
baryons in the Universe. Since many (most?) baryons 
may be "dark" , such approaches can best set lower 
bounds to the present ratio of baryons-to-photons. One 
such estimate 22 suggests a very weak lower bound on 
77 of: 7710 > 0.25, entirel y c onsistent with the BBN es- 
timates above. Others |2^ have used more subjective 
(although cautious) estimates of the uncertainties, find- 
ing a much higher lower bound to the global budget of 
baryons: 7710 > 1.5, which is still consistent with the 
"low-77" range identified using the high-D results. 

A possible challenge to the "low-7;" case comes 
from an analysis |24| which employed observational con- 
straints on the Hubble parameter, the age of the Uni- 
verse, the structure-formation "shape" parameter, and 
the X-ray cluster gas fraction to provide non-BBN con- 
straints on the present density of baryons, finding that 
7710 > 5 may be favored over 7^10 < 2. Even so, a signif- 
icant low-77, high-D range still survives. 

5 Cosmology Constrains Particle 
Physics 

Limits on particle physics beyond its standard model 
are mostly sensitive to the bounds imposed on the '^He 
abundance. As described earlier, the "^He abundance 
is predominantly determined by the neutron-to-proton 
ratio just prior to nucleosynthesis. This ratio is deter- 
mined by the competition between the weak interac- 
tion rates and the universal expansion rate. The latter 
can be modified from its standard model prediction by 
the presence of "new" particles beyond those known or 
expected on the basis of the standard model of parti- 
cle physics. For example, additional neutrino "flavors" 
{AN„ > 0), or other new particles, would increase the 
total energy density of the Universe, thus increasing the 
expansion rate (see equations 4 & 5), leaving more neu- 
trons to form more ''He. For ANi, sufficiently small, the 
predicted primordial helium abundance scales nearly 
linearly with AN^: AY « 0.013AiV^. As a result, an 
upper bound to Yp coupled with a lower bound to 77 
(since Yp increases with increasing baryon abundance) 
will lead to an upper bound to AiV^ and a constraint 
on particle physics Q. 

The constraints on N^, = 3-1- AA''^, as a function of 
the baryon-to-photon ratio 77 is shown in Figure 5. For 
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Figure 5: The number of equivalent massless neutri- 
nos Ni, = 3 + AN^ (see eq. 4) as a function of rj. The 
region aUowed by an assumed primordial ^He mass frac- 
tion in the range 0.228 - 0.248 lies between the two solid 
curves. The vertical bands bounded by the dashed lines 
show generous bounds on rj from "high" and "low" deu- 
terium. 

"low-D/high-77" there is little room for any "extra" par- 
ticles: AA''^ < 0.3. This would eliminate a new neutrino 
flavor {ANi, = 1) or a new scalar particle (AiV^ = 4/7) 
provided they were massless or light (m ^ 1 MeV), and 
interacted at least as strongly as the "ordinary" neutri- 
nos. In contrast, much weaker constraints are found for 
"high-D/low-r;" where N^, as large as 5 {ANi, < 2) may 
be allowed (see Fig. 5). 

6 SUMMARY 

The standard, hot, big bang cosmological model 
is simple (assuming isotropy, homogeneity, Newto- 



nian/Einsteinian gravity, standard particle physics, 
etc.) and, likely, simplistic. In broad brush it offers 
a remarkably successful framework for understanding 
observations of the present and recent Universe. It may 
seem hubris to expect that this model could provide 
a realistic description of the Universe during the first 
20 minutes or so in its evolution when the tempera- 
ture and density were enormously larger than today. 
According to the standard model, during these first 20 
minutes the entire Universe was a nuclear reactor, turn- 
ing neutrons and protons into the light nuclides. This 
prediction presents the opportunity to use observations 
here and now to test the theory there and then. As 
described in this article, SBBN predicts observable pri- 
mordial abundances for just four light nuclides D, ^He, 
^He, and ^Li, as a function of only one adjustable pa- 
rameter, rj, the nucleon-to-photon ratio, which is a mea- 
sure of the universal baryon abundance. Even though 
it is currently difficult to use the extant observational 
data to bound the primordial abundance of '^He, SBBN 
is still overconstrained, yielding three predicted abun- 
dances for one free parameter. Furthermore, the baryon 
density inferred from SBBN must be consistent with 
that derived from observations of the present Universe. 
Given the many possibilities that SBBN could be falsi- 
fied by the empirical data, it is a remarkable success of 
the standard model that there is consistency between 
theory and observations provided that there are a few 
billion photons (most of them in the 2.7 K cosmic back- 
ground radiation) for every neutron or proton (nuclcon) 
in the Universe. 

This success establishes primordial nucleosynthesis 
as one of the main pillars of our standard model of cos- 
mology, providing the only probe of the physical Uni- 
verse during its very early evolution. Alternate theories 
of gravity and/or particle physics must now be tested 
against the impressive success of SBBN. 
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